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Abstract 

Low-temperature dynamics of flux lines in high temperature, type II, super- 
conductors in the presence of correlated disorder in the form of columnar de- 
fects is discussed. The effect of tilting the applied magnetic field with respect 
to the column's directions is considered, using the non-Hermitian quantum 
mechanics technique evaluated by Hatano and Nelson It is shown that 
the critical current, as well as the vortex transport properties below this cur- 
rent, may be determined by "surface excitations", i.e., by the roughness of 
the flux line near the edges of the sample, which dominated the bulk jumps. 
Phase space considerations determine the critical thickness of the sample, 
below which the tilt induced surface transport exceeds the bulk mechanism. 
This critical length, which depend on the tilt angle as well as the directions 
of the perpendicular field and the suppercurrent, diverge at the derealization 
transition. 
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Flux lines response functions in cuperate high-temperature superconductors have at- 
tracted considerable interest in recent years 0. In order to avoid dissipation of energy as 
a result of flux lines motion driven by the superconducting current, these lines should be 
pinned by crystal impurities ||. It turned out |4]] that the pinning is much stronger (es- 
pecially when many vortex interaction is taken into account) when these impurities are in 
form of correlated disorder (such as twin boundaries f| or columnar defects ||) compared 
to the case of point disorder, resulting from vacancies of Oxygen atoms [|7| . However, Nelson 
and Vinokur Q have pointed out that the correlated defects pinning become less effective in 
cases where the direction of the external magnetic field is tilted with respect to the defects, 
which we take to be along the z direction. At some critical tilt, for which the energy per 
unit length of the defect is less than the energy associated with the perpendicular field, a 
pinning-depinning phase transition occur and the flux lines delocalized. 

Critical bulk current and vortex dynamics below this current for flux lines in the presence 
of columnar defects have been considered in [|J. The authors, using the mapping of flux lines 
in d+ 1 dimensional superconductor to the world lines of bosons in a d— dimensional quantum 
system, identified the phase space diagram of the system which contain a high temperature 
"superfluid" and low-temperature "bose glass" phases, as well as Mott insulator at the 
matching field, = n P i n (po, for which there is one flux line per defect. At low temperature, 
this matching field separates the "dilute" region of the bose glass phase, for which the vortex 
lines are pinned individually by the defects, (i.e., a ~ (4> /By/ 2 , the Abrikosov lattice 
constant, is much bigger than d, the typical distance between two columnar defects) from 
the high density region, for which interactions are important in determining the localization 
length and transport properties of the flux lines. 

In the low field region, the vortices are localized by the interaction with the correlated 
defects. Each defect is the analog of a 2D potential well which we ??? take as a cylindrical 
square well such as V(r) = —Uq for r < bo and V(r) = for r > bo- The temperature of the 
superconductor, in turn, corresponds to the Plank constant h of the quantum boson problem. 
For the dilute vortex arrays, where the pinning energy is larger then the interaction energy, 
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there are two regimes. For low temperatures (T << T* , T* = \/Uo^ibo) the localization 
length l± is on the radius of the defect; i.e., of order 6 ; so that each flux line is localized 
by one defect. On the other hand, for T » T*, the localization length of one defect grows 
exponentially with T 2 , and the flux line is then localized be several defects, forming an 
effective d dimensional potential well in the corresponding boson system. 

The response of the flux line to superconducting current in the plane perpendicular to the 
vortex direction J _L B translates itself into the response of the boson system to an applied 
electric field. For vortices oriented in the z direction, the Lorentz force per unit length of 
the vortex is given by 

f L = x J (1) 

c 

Which is the analog of a boson with charge 0o interacting with the electric field E = 
±z x J. 

c 

Above the critical current J c , the vortices are no more localized ant there is no supercon- 
ductivity (in the sense of dissipation free current) anymore. Below this critical current, the 
mechanism for flux flow is tunneling via thermally activated "half loops" (or, in the boson 
dynamics, tunneling into the conduction bend). For currents smaller than ji, the half loop 
transverse displacement exceeds the mean distance between occupied pinning sites, and for 
thick sample the flux lines move via the nucleation of superkinks, the analog of the Mott 
variable range hopping (VRH) in doped semiconductors 0. 

The depinning of the flux line as a result of external field tilt has been carefully inves- 
tigated by Hatano and Nelson |IJ. The Hemiltonian of the corresponding boson problem is 
no longer Harmitian; the kinetic term of the Hamiltonian p 2 /(2ei) (p = — zTV) is subject 
to imaginary gauge transformation and takes the form (p + ih) 2 /(2ei), where h is related 
to the perpendicular magnetic field Hj_ via h = </> H_|_/(47r). As a result, there are two 
solutions for each localized (real energy spectrum) state, corresponding to the eigenstates of 
the Hamiltonian and its complex conjugate. These solutions, termed ipR and ipL, correspond 
to the right or left "tilting" of the localized solution of the untilted Hamiltonian, i.e., 
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^L(r)=exp(±h.r)^(r). (2) 

The probability distribution to find the flux line at the point r at a distance r from the 
surface of the sample is given by 

P(t,t) = Z- 1 < ip f \ exp(-(L -t)H/T) |r >< r| exp(-r#/T)|^ > (3) 
= Z- 1 J2< ¥\m L >< m R \r >< r\n L >< Ur^ > e -(rE m +(L- T )E n )/T ^ 

m,n 

such that, as L — > oo, the probability distribution of the flux line at the surface is propor- 
tional to < r\gs, L, R >= ^ R s{v) where \gs L>R > are the left and the right ground state, 
respectively. Deep in the bulk, the distribution is given by P(r, L/2) = ip g s R ipg S = ipuntutedi 
i.e., in the localized regime, the flux line changes its shape near the surface, while remains 
uneffected in the bulk (see Fig. 1). Typically, the "surface roughness associated with the 
tilt extends into the bulk up to characteristic distance r*, which diverges as the tilting angle 
approaches the critical angle, for which the flux line delocalizes and the current response 
become linear. 

In this communication, we study the the effect of the tilt on the 

flux lines response to superconducting currents in the regime where the tilting angle is 
smaller then critical, i.e., in the bose glass phase where the vortex lines are localized. We 
assume that the thickness of the sample is large enough, such that it is much larger than the 
dimension of the optimal excitation along the z exis, Moreover, we address only the dilute 
limit, for which the transverse (xy) displacement of these excitations is less than a , so that 
the interaction is taken into account by filling up the localized states in order of increasing 
energy up to the chemical potential /i. 

Let us consider first the critical currents. This current is determined by the binding free 
energy U(T) as well as the localization length l± = 1/k. Modeling the defect as a square 
well in the boson system, k is related to U(T) by k = ^2U{T)e\/T . Of these two, l± is 
changed as the magnetic field is tilted. Near the surface of the sample, the localization 
length should be l±(h,9) = 1/(k — h\cos(9)\), where 9 is the angle between the Lorentz 
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force fi (perpendicular to the supercurrent J, which we take to be in the xy plane) and 
the tilting field h. The absolute value is needed for the case of \9\ > 7iover2, for which the 
critical current is dominated by the "tail" of the flux line on the other end of the sample, as 
shown in Fig. 1. Thus, the critical current will take the form, 

fflM) = g^Wf. (5) 

00 

This critical current is determined by the surface ends of the vortex, for which the effect of 
the tilt is maximal. However, the "creep" of the vortex in the direction of the tilting field is 
limited by the effect of "image vortices" which should be introduced in order to satisfy the 
boundary condition on the surface |4[]. Thus, for very small perpendicular magnetic field, 
where the surface roughness extension r* is less then the London length, the tilt have no 
effect on the vortex pinning and there is no change in the response functions of the flux 
system. 

For currents less then critical, the thermally assisted flux flow (TAFF) theory of the 
vortex transport gives the resistivity p = £ / J as 

p = poe 5F ' T (6) 

where SF is the energy barrier for flux line jumps. Our basic observation is that deep in the 
bulk there is no influence of the tilt, so that the energy barriers for nucleating half loops or 
double kinks are the same. The physical reason for it is that, although the perpendicular 
field decreases the energy barrier for one side of the kink / loop, it increases the energy 
needed for the other side. The main effect of the tilt comes from surface kinks / loops, for 
which the energy barrier really decreases. Although the resulting energy barrier SF is smaller 
then the bulk one, so that the "resistivity" associated with it is exponentially smaller, one 
should take into account the phase space prefactor of these two mechanisms - the number of 
surface kinks available is determined by the surface roughness, i.e., by r*, while the number 
of bulk kinks is of order (L — t*)/Z, where Z is the distance for which the half loop / kink 
extends along the relevant defect. Thus the nature of the current response is determined 



by the thickness of the sample - for L > L c (h), the bulk excitations will dominate and the 
response is tilt independent, while for L < L c surface excitations become important and the 
voltage drop will be tilt dependent. As h — > h c (where h c is the critical field above which 
the flux lines delocalize) L c diverges, so that near the depinning transition the resistivity of 
the system goes continuesly to zero. 

In order to estimate the relevant quantities in the tilted case, we use the expression for 
the free energy of the surface excitations in the presence of the tilt. Consider now surface 
excitation of the flux line with line tension e\ which extends for a distance z along the pin 
and has perpendicular extent r. The free energy of such jump is given by 
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5F = — + U z - f L rz - hr cos(#) (7) 

for the "half loop" surface excitations, if the jump is due to the nucleation of superkinks, 
one should take into account the energy difference between different rods at distance r. This, 
in turn, is determined by the density of states at the chemical potential [Q. and the 
free energy is, 
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5F= 6 -^- + — ^— - f L rz - hr cos(0) (8) 
z g(i2)r 2 

The resulting saddle point free energies are 

6F* = (E k - h\ cos(0)|d)(Ji/J) halfloops (9) 

5F* = (E k - h\ cos(6)\d)(J /J) 1/3 superkinks (10) 

where E k = y/eiUod, J\ = cUo/(4>od) and Jo = cj (4>og(fi)d 3 ), for d the average spacing 
between unoccupied pins. 

Let us estimate now the phase space for such surface excitations, i.e., the width of the 
region in which the roughness takes place. Using Eq. (3) one finds that the crossover 
between the surface (P(r) ~ ip R,L (r)) and the bulk, for which P(r) is the same for the tilted 
and the untilted situation, is determined by the quantity 
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K(r) = ]T < ™ L \r > ex P (-rE m /T) (11) 

m 

where r is the distance from the surface. The transition to the surface behavior takes place 
when Y(r) becomes r independent, and thus absorbed into the normalization factor for 
P(r). Typically, this happens when (10) is not dominated by the r dependent exponential 
factor, since then the summation over m is determined by the delocalized states, yielding 
an r independent result. Thus, for the case of half loop tunneling, the width of the surface 
roughness will be r*(h,6) « T/E*(h,9) where Tk(E*) = h. This gives us the estimate 

* H/M 

loops - _ h 2y ^ ) 

For the superkinks tunneling, the energy i? m may be given by l/(g(/i)r^J, so that the 
energy exponent become negligible as r << r*, where [|TJ, 

T = JK-^W (13) 

The phase space of the surface excitations is given by the width of the surface region 
divided by the "width" of the typical excitation z* ur f ace . Using the above expressions for 
the free energy of the kinks/loops, one finds that 

* (Jhe) c(y/^Z7o-%Os(fl)|) 

^surface,loops\°i"'i v )~ ji \ L ^J 

J(p 

and 

{J ) h ) 6) = c{Ek -^ m (15) 
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sur face, kinks \" i ' ~i ~ i ji ^ 

The resulting resistivity in thick samples will be determined by adding in parallel the 
t* I z* "surface resistors" with p = e SFsur f ace ^ T to the system of (L — r*) / Z£ ulk (J) "bulk 
resistors" with Z£ ulk = z* ur j- ace (h = 0). While the surface roughness does not depend on 
the angle between the current and the transverse magnetic field, the width of the jump, 
as well as the free energy barrier, do depend on it. It turns out that the resistivity in the 
"perpendicular" direction (ft _L h ± z) is independent of the tilt. For other directions of 
the superconducting current, there will be a crossover length L c below which the surface 



loops dominate the jumps. For any tilt less then critical, the surface roughness is finite, 
so that as L — > oo bulk excitations are clearly the preferred jumping mechanism, but as 
h — > h c , the width of the surface roughness become comparable with the sample thickness 
for each finite sample and one sees a crossover to surface excitation dominated transport. 
The critical length is related to the parameters above as 

L c ( J, h, 9) = r *(_^l_ exp(^ fc - 5F: urface ) + 1). (16) 

surface 

There are two reasons for the divergence of L c as h — > h c ; one is the divergence of 
r*, the other is the fact that z* ur j- ace — > 0, yielding infinite phase space for the surface 
excitations. However, there is a limitation on the minimal width z of the jumps; as z toX, 
the London length, self interaction of the flux line lock the kink/loop, so that A sets the 
minimal excitation extent along the z axis. For the region in parameter space for which 
z* » A, the critical length will grow like (h\cos(9)\ — /i c )~ 2 for loop transport, and as 
(h\ cos(0)| — h c )~ 3 for kinks. On the other hand, as one approach the critical tilt, the region 

(^-^08(0)1) <^ (17) 

is entered, in which the excitations width could not shrink anymore. In that case the critical 
length diverges like (h\ cos(6>)| — /i c ) -1 for loops, and as (h\ cos(0)| — h c )~ 2 for kinks. 
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Figure Captions 

Figure l:Flux lines locaqlized by columnar defects in the presence of perpendicular field 
Hj_. The surface roughness extends distance r into the bulk, and the Lorentz force fi, is at 
angle 9 to the tilting field. For \9\ > ir/2, the contribution of the surface roughness to the 
vortex transport comes from the "tail" at the lower end of the sample in the figure. 
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